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■ Abstract. In this paper we study how to accelerate the convergence of the ra- 
04 ' tios (xn) of generalized Fibonacci sequences. In particular, we provide recurrent 

(— j , formulas in order to generate subsequences (x Bn ) for every linear recurrent se- 

■ quence (g n ) of order 2. Using these formulas we prove that some approximation 

methods, as secant, Newton, Halley and Householder methods, can generate sub- 

. sequences of Moreover, interesting properties on Fibonacci numbers arise as 

an application. Finally, we apply all the results to the convergents of a particular 
continued fraction which represents quadratic irrationalities. 

H 

> ( 

1. Introduction 

The problem of finding sequences of approximations for a solution of a certain 
equation f(t) = is really intriguing. There are various and famous methods to 
I generate such sequences in literature, as Newton, secant and Halley methods. But 

C^- ■ another interesting problem is how to increase the rate of convergence for the dif- 

, ferent methods. One of the answer consists in selecting some subsequences of the 

\ starting approximations sequences which accelerate the convergence process. In the 

particular case when f(t) corresponds to a second degree polynomial, we may think 
I about it as a characteristic polynomial of a generalized Fibonacci sequence. The ra- 

' tio of its consecutive terms generates a new sequence which converges to the root of 

larger modulus of f(t), when it exists and it is real. Many authors have given meth- 
ods to accelerate the rate of convergence of this sequence. McCabe and Phillips [5] 
have studied how some methods of approximations, like Newton and secant meth- 
5_i , ods, provide an acceleration of such sequence for suitable initial values. Gill and 



Miller [2j have found similar results for the Newton method, only for the ratios of 
consecutive Fibonacci numbers. Many other works have been developed about this 
argument. In [4] and [6j, e.g., the results of McCabe and Phillips [5] have been gener- 
alized for ratios of non-consecutive elements of generalized Fibonacci sequences. In 
[7j these results have been extended to other approximation methods, like the Halley 
method. The aim of this paper is to give new interesting points of view to prove 
in an easier way some known results about accelerations of generalized Fibonacci 
sequences, providing recurrent formulas in order to generate these subsequences. 
The accelerations related to approximation methods become particular cases of our 
work. In the first section we will expose our approach, showing in the second sec- 
tion the relationship with known approximation methods as Newton, secant, Halley 
and Householder methods. Moreover our results, in addition to a straightforward 
proof of some identities about Fibonacci numbers, can be applied to convergents of 
a particular continued fraction representing quadratic irrationalities. We will prove 
in the last section that particular subsequences of these convergents correspond to 
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Newton, Halley and secant approximations for the root of larger modulus of a second 
degree polynomial, when it exists and it is real. 

2. Accelerations of linear recurrent sequences of the second order 

Definition 2.1. We define, over an integral domain with unit R, the set S(R) of 
sequences a = (o n )n^o> an< ^ ^ e se ^ VV(i?) of linear recurrent sequences. Further- 
more, we indicate with a = (a n )^= = W(a, /3,p,q) the linear recurrent sequence of 
order 2 with characteristic polynomial t 2 — pt + q and initial conditions a and j3, 
i.e., 

oq = a 

ai = (3 

a n = pa n _i - <?a n _ 2 Vn > 2 . 

If we consider the sequence U = W(0, l,p,q), it is well-known that the sequence 
x = (x n )^= 2 defined by 

x n = 7^-, Vn > 2 (2.1) 
Un-l 

converges to the root of larger modulus (when it is real and it exists) of the char- 
acteristic polynomial of U. Thus, every subsequence {x gn )^^ (for some sequence 
(g n )n=o) can b e considered as a convergence acceleration of the starting sequence 
x. Cerruti and Vaccarino [l] have deeply studied the applications of the companion 
matrix of linear recurrent sequences. Here we consider some of these results which 
are useful to find the explicit formulas for the acceleration of the sequences U and 
x. We define 

' o v 

-q p ; 

as the companion matrix of a linear recurrent sequence with characteristic polyno- 
mial t 2 — pt + q. We recall that some authors may use the transpose of M as the 
companion matrix. 

Lemma 2.2. The sequences U = W(0,l,p,q) and T = W(l,0,p,q) satisfy the 
following relations 

~ qUn , Vn>0. (2.2) 

: T n + pU n 

\l" - I '" Un I = ( Tn Un 1 (ID 

~ \T n+1 U n+1 J \-qU n T n +pU n )' {Z -°> 

where M is the companion matrix of U and T. 

Proof. Using the Binet formula and initial conditions, we easily check that 

a n - /3 n afi n - pa 11 

U n = 3-, T n = — , Vn > 0, 

a — p a — p 

where a and f3 are the zeros of the characteristic polynomial t 2 — pt + q. Thus, 
remembering that a/3 = q, we immediately obtain (|2.2f) . In order to prove (12 .3p . we 
observe that the characteristic polynomial of U and T annihilates M, i.e., we have 

1 0\ , , / 1\ , ,,, . ^n— 1 m „ /rn _2 



M 




= ' 1 J ' M = [ -q p ' ' ^ = vM ' n ~ ~ qMU ~ ' Vn - 2 ' 
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So, using (j2.2|) . we clearly find that the entries of M n in the first column are the 
terms T n , T n+ \ of T, and in the second column are the terms U n , U n +\ of U. □ 

Dealing with definition (|2.ip of x, the ratio of two consecutive terms of any sequence 
a € W(R) has a closed formula, as we show in the following 

Theorem 2.3. Let us consider any sequence a = W(ao,ai,p,q) and the sequence 

V = {Vn)n=2> with Vn = — —, then 

a n a x x n ~a q 

Vn = = ; Vn > 2 , (2.4) 

a n -\ a x n + ai - a p 

where the n-th term of the sequence x defined by 112.1]) . 

Proof. By direct calculation we have 

a n aiU n + a T n a\U n + a T n 

Vn 



a n -i aiU n -i + a T n _i a\\J n -\ + a (pU n -i - qU n - 2 - pU n -i) 

Un 

aiU n - a qU n -i ai \J n _ x a ° q 



aiU n ^i + a U n - a pU n -i U n 

a — h a\ - aop 

Un-l 

aix n - a q 



aox n + ai - aop 

□ 

Now, as a very important consequence of the previous theorem, we can evaluate the 
terms of the sequence x shifted by m positions 

Corollary 2.4. If x is the sequence defined by \2.1\) . then for all m > 2 — n 

X n + X m +i p 

Proof. Let us consider a = W(U m , U m +i,p,q), that is a n = U n + m for all n > 0. So 
the sequence y, introduced in the previous theorem, becomes y = (y n = x n+m )^^ 2 ^ 
and, from (|2.4p . we obtain 

U m +i%n U m q x m -\-ix n q 

Xn+m = 77 —f7 77 = ; • (^-5) 

U m x n + U m+ i - U m p x n + x m+1 -p 

□ 

The basis sequences U and T play a central role in the following results. The acceler- 
ations for such sequences yield explicit formulas for the acceleration of the convergent 
sequence x = (a? n )l^. Moreover, some interesting properties on Fibonacci numbers 
arise as an application. We recall without proof an important theorem about the 
generalized Fibonacci sequences. 

Theorem 2.5 (pQ, Th 3.1 and Cor. 3.3). Denote by V = W(2,p,p,q) the Lucas 
sequence with parameters p and q, then 

W mn (0, l,p, q) = W m (0, l,p, q) ■ W n (0, 1, V m ,q m ), 
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W mn {l, 0,p, q) = W n (l, 0, V m , q m ) + W m (l, 0,p, q)W n (0, 1, V m ,q m ). 

Theorem 2.6. Given the sequences U = W(0, l,p, q), T = W(l, 0,p, q), and x as 
defined in equation \2. 1\) . if we choose a sequence g = W(i,j, s, t), for some i, j > 2, 
then 



U 9n = a 2 U^_ 1 T^l + hT^U^l + ( ai b 2 + a 2 6 4 )C/W_ 1 C/j, 



-*) 

9n-2 



for every n>2, where = W(0, 1, V s , q s ), = W(l, 0, V s , q s ) and 

M s = (a l a 2 \ M -t = (h b 2 
\a 3 a 4 J \b 3 6 4 

Moreover, if we choose Xi and Xj, for some i, j > 2, as the initial steps of the 
acceleration of x, then 

gWrff., + qH 2 x g - n % - q( ai b 2 + a^^t^t-l 

Xgn n2 nn J') nh J~t) ,( nh , _ u W («) _(-*) 
q z - gaix 9n _ 1 - qo\x 9n _ 2 + [a\b\ + a^jXg^Xg^ 

, . [/« 
where x (s> = —q — rr • 

Proof. We know as a consequence of Theorem 12.51 that 



Vn > 2 , (2.6) 



M S9n ~ 



!\/r-*9n-2 _ (Tgn-l + b\Ug n % b 2 U 9n % 



Tg n ^j?n \ _ j^9n — Jy[Sgn-l—tg„-2 — ]\/[ s 9n - 1 — *9n-2 



So from 

we easily find the explicit expressions of ?7 9n and T 9n . Finally, to complete the proof, 

we only need to divide both U 9n and T 9n by the product Tg^^Tg^ and consider 
their resulting ratio multiplied by — q 

T T (s) rr(-t) rr(s) rr(-t) 

n- a2 ^T +62 ^ +(aib2 + a2fe4) T M T (-t) 

_ ^ffra _ _ £gn-l ± gn-2 £gn-l ± 9 n-2 

1 , f™- 1 . , u !?n-2 . / , . , v u g n -i u g n -2 

I g n -i I g n -2 1 g n -i I g n -2 

rearranging the terms we obtain the equality (|2.6p . D 

Corollary 2.7. WTien 3 = W(2, 3, 1,-1), i.e., £r n _ 3 = F n , Vn > 3, where F = 
W(0, 1, 1, —1) is the Fibonacci sequence, then for all n > 5 

UF n = Up^Tp^ + T Fn ^ U Fn „ 2 + pU Fn _ x U Fn _ 2 = 

= ~Q (f/F n _!^F n _2-l + Up^-lUp^) +pU Fn _ 1 U Fn _ 2 , 

qx Fn _, +qx Fn _ 2 -px Fn _ x x Fn _ 2 



q-x Fn _ x x Fn _ 2 
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Proof. The identities immediately follow from Theorem 12.61 observing that in this 
case we have s = 1, t = —1, consequently M s = M = M~*, = U = £A - *) and 

rp{s) — J< — ^(-t). □ 

Remark 2.8. When U is the Fibonacci sequence F, the previous Corollary provides 
a direct proof of the interesting formula 

F Fn = Fjr^Fp^! + Ff^-iFf^ + Ff^Ff^ , 

for all n > 5, avoiding many algebraic manipulations based on the basic properties 
of Fibonacci numbers. This formula is also true for n = 3 and n = 4, as a little 
calculation can show. 

The acceleration shown in the previous corollary is the same found by McCabe 
and Phillips in [5] with the secant method applied to the sequence (x n )^= 2 shifted 
by one step. 
If we consider 



M 2n = M n M n 



Tl-qVl 2T n U n +pU% 
-q{2T n U n +pUl) Tl-qUl+p{2T n U n +pUl)) ' 

we find U 2n = 2T n U n + pU 2 , T 2n = T% - qU 2 n , and 



X2n = -It^t- = — : 2 — • <- 2 - 7 ) 



Ihn _ 2qx n - p x\ 

T 2n q- x\ 

If we start from x 2 , a repeated use of equation (12. 7p yields to the subsequence 
x 2 , X4, xs, x 2 n, corresponding to the subsequence obtained in [5], using the 
Newton method on the sequence (i n )^™ 2 shifted by one step. 

Proposition 2.9. Using the notation of Theorem \2.6\ for g = W(h,h + k,2,l), 
i.e., g n = kn + h Vn > 0, then for all n > 2 

U * = i U U l-i U 9 n -* + *T 9n -iU 9n _,T gn _ 2 +p Ul_ 1 T 9n _ 2 - U 9n _ 2 Tl^ 



and 



Tg n — q 9n _ 2 (^g„_i^9n-2 +pTg n l Ug n _ 2 QT 9n _ 2 U gn _ 1 + 2qT gn _ 1 U gn _ 1 U gn _ 

X ln-l X 9n-2 + 2qXg n _ 1 ~ P% g „_ 1 ~ Q X g n - 2 



1 P X 9n-2 X< g n -1 2,X 9n-l X gn-2 

Proof. We need to evaluate 

M 9n = M 2g,l - 1 M~ gn - 2 . 



Considering that 

M~ 

we have 



q n \q J ~ q n\ qU n T n 



= — ~ « U l-i 2T 9n-.U 9n ^ +pU 2 9n _ 1 \ (T gn _ 2 + pU gn _ 2 -U gn _ 2 
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and 

U *> = ^ {l U L^ 9 n-, + ZT^U^T^ + P Ul_ 1 T gn _ 2 - U^T 2 ^ 

T 9n = ^{ T Li T ^-> +PTL 1 U gn ^-qT gn _ 2 Ul_ 1 +2qT gn _ 1 U gn _ 1 U gn _ 2 
Finally, dividing both U 9n and T gn by Tg n i T gn _ 2 , their ratio becomes 

^9n-l \ Ugn-2 | 2 ^ 9n ~ 1 | p( ^ gn ~ 1 \ ^3n-2 
Ug n V^9n-1/ F g n -2 Tg n _ 1 \ F g n -l J Tg n _ 2 



1 +P7^ 9 Tfi + 2 <? 



± 9n-2 \ ± g,i-i/ ± gn-\ ± g n -2 

Ug n 

from which, remembering that x gn = —q——, with simple calculations we obtain the 

thesis. □ 

Remark 2.10. Applying the last Proposition to the Fibonacci numbers, we get an- 
other example of a beautiful identity, easily proved without hard calculations. In 
fact, when we consider g n = kn, we have for all n>2 

F kn = {—^) {— F k{n-l) F k{n-2) + ^ F k(n~l) F k(n-l)-l F k{n-2)-l 

+ F k{n-l) F k(n-2)-l ~ F k{n-2) F k{n-l)-l) ' 

which in particular , for k = 1, becomes 

F n = ( — I)" {— F n-l F n-2 + 2 F n-l F n-2 F n-3 + F n-l F n-3 ~ F n-2) ■ 

3. Approximation methods 

In the previous section we have already seen how some accelerations are connected 
to some approximation methods, like Newton and secant methods. In this section, as 
particular cases of our work, we exactly retrieve the accelerations provided by these 
methods and studied in the papers [2j and [5] . Furthermore, using our approach, we 
can show similar results for different approximation methods, like the Halley method. 
Finally, we will examine the acceleration provided by Householder method, which 
is a generalization of Newton and Halley methods. First of all, we recall that the 
secant method, applied to the equation f{t) = 0, provides rational approximations of 
a root, starting from two approximations xq and x%, through the recurrence relation 

_ f{x n -l)(x n -i - X n - 2 ) 

Xfi — X n — 1 . . . . , VTJ Z, 

f{x n -l) ~ f{x n -2) 

which becomes, when fit) = at 2 — bt — c 

ax n -ix n -2 - c 



x , 



ax n _\ + ax n -2 — b 



Theorem 3.1. Let us consider the Fibonacci sequence F and the sequences U = 
W(0, l,p, q), T = W(l, 0,p, q), and x = (x n )'^5 2 , defined by equation \2. The 
subsequence {xF n+2 +i)n=o corresponds to the approximations of the root of larger 
modulus of t 2 — pt + q provided by secant method. 
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Proof. We take the sequence g defined by g n = F n+ 2 + 1 for every n > 0. Since 
g n — 1 = -Fn+2, we have the recurrent formula g n = g n -i + <?n-2 — 1. Let M be the 
companion matrix of the sequence U. From the matrix product 

M 9n = M 9n - 1 M 9n - 2 M~ 1 

we immediately observe that 

1 p 
tj — J j tj T T T — U T A- U T 4- —T T 

U 9n — U g n -1 U 9n-2 J-9n-\ ± 9n-21 ± 9n ~ U 9n-2 ± Sn-1 ' U ffn- 1 1 9n-2 ' * 9n- 1 ± g-n-2 ' 



Thus 



U T T 

U 9n - t g n _l- t 3„_2 



T 9n Ug n _ 2 ^ Ug^ p 



and 



T T 



Ug n 



y n = y n _i H ^ — — ^ — - — , Vn > 1 , (3.2) 



X 3n-2 + ^Sn-l P 

clearly proving the thesis. □ 

The Newton method provides rational approximations for a solution of the equa- 
tion f(t) = 0. When f(t) = at 2 — bt — c, given an initial approximation uq, we recall 
that the Newton approximations sequence satisfies the recurrence relation 

ayl_ x - by n -i - c ay 2 n _ x + c 

Vn = Vn-1 t, r = o 7, V n > 1 . (3.1) 

2oy n _i - b 2ay n -i - b 

On the other hand, the Halley method generates approximations for a solution of 
t 2 — pt + q = through the following recurrence relation 

— PVn-l + q){p- 2j/n- l) 
3 2/n-l - 3 Pyn-l +P 2 - 9 

for a given initial step yo- I n the next theorem we will show how particular sub- 
sequences of the sequence x, defined by equation (|2.ip . generate the Newton and 
Halley approximations to the root of larger modulus of t 2 — pt + q . 

Theorem 3.2. Given the sequences U = W(0, l,p, q), T = W(l,0,p, q), and x, 
defined by equation (2.1}) . then the recurrence relations 

x 2 n — q 

X2n-1 = Tj , (3-3) 

2x n -p 

, (»n_l -pXn-l +q){p- 2x„_i) 

X 3 n-2 = + 7T^ 2 q 1 2 ' ( 3 " 4 ) 

3<_i - 3px n _i - q 
respectively generates the subsequences of x producing Newton and Halley approxi- 
mations for the root of larger modulus of t 2 — pt + q . 



Proof. The proof of relation f|3.3|) is a direct consequence of Corollary [27 



X2n—1 — X n -\- n — i — 

2x n -p 
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Starting from x 2 , a repeated use of this equation yields the subsequence 

(X 2 ,X 3 ,X 5 , . . . ,X 2 n + l, ■ ■ •) , 

which consists of the Newton approximations for the root of larger modulus of 
t 2 — pt + q . Relation <\3.3h and Corollary 12.41 also provide the proof of recurrence 
(|3.4p for Halley approximations 

X2n-\X n ~ q 

— X n+ 2n-2 



X n + X 2n -1 ~ P 

2x n - p X " ~ q x 3 n - 3qx n + pq 



2 

x n — q 



x 2 n -q 3j;2 _ 3 pXn + p 2 _ q 

Xn + t, V 

2x n -p 



{xi_i - px n -i +q){p- 2x 



n-ll 



3xt l _ 1 - 3px n _i +p 2 - q 



□ 



To complete our overview on approximation methods and generalize the results of 
previous theorem, we consider the Householder method [3] , which provides approxi- 
mations for a solution of a nonlinear equation fit) = 0. The Householder recurrence 
relation is 

J (i//) |J - 1> (^) ., . _ 

+ i (!//)(«%„) ' V " a ° (3 ' 5) 

for a given initial approximation y®, where /W is the d-th derivative of /. The 
Newton and Halley method clearly are particular cases of this method, for d = 1 
and d = 2 respectively. Indeed, when fit) = t 2 — pt + q from (|3.5p . when d = 1 or 
d = 2, it is easy to retrieve the formulas (|3.1|) and 



Theorem 3.3. With the scone hypotheses of Thcoreui \3. S\ , if y n — for some 
k > 2, then y n +i = XM+iU-dj where y n +i is obtained from relation (3. 5\) when 
f(t) = t 2 -pt + q. 

Proof. Let a±, c*2 be the roots of t 2 — pt + q, we can write 

/ 1 \ (d) _ ( (t - ai )- x - (t - a 2 )-!\ {d) _ {-l) d d\ ((* - a 2 ) d +! - (t - ai ) d+1 ) 



yt 2 -pt + qj y ai-a 2 J {a 1 -a 2 )(t-a 1 ) d + 1 (t-a 2 ) d + 1 

Thus, relation (|3.5p becomes 



and with some algebraic manipulations we obtain 

ai(y n ~ a 2 ) d+1 - a 2 (y n - 



Vn+1 



d+l 



ACCELERATIONS OF GENERALIZED FIBONACCI SEQUENCES 



I) 



Uk oq — oQ 
If we set y n = x k = j- = -ett uZT, then 



Vn+l 



(u k-l\ d+1 / fc-1 fc\ d+1 

al - a 2 a^ \ / a\a 2 —"21 
~k=i kZT - «2 fcTi 
Q=l - «2 / \f*l ~°2 / 

(fc k-l\ d+1 / fc-1 k\ d+1 

af - a%a\ \ I a\a 2 -02) 

k-l 1 ~~ [ k-l k-l J 
a x -a 2 J \a x - a 2 J 

a (d+i)k-d {ai _ a2)d+1 _ a (d+i)k-d( ai _ a2) d+i 



(d+l)fc— d— 1 / njii (d+l)fc— Nj.i 

1 (ai — 02) — «2 ( a i ~~ «2j 

^(d+ljfc-d • 

□ 



U(d+l)k-d 
U(d+l)k-d-l 



4. Applications to continued fractions 

In order to show some interesting applications of the exposed results, we study 
the acceleration of a sequence of rationals which come from the sequence of con- 
vergents of a certain continued fraction. The continued fraction, which we will 
introduce, provides a periodic representation of period 2 for every quadratic irra- 
tionality (except for the square roots). Surely its convergents do not generate the 
best approximations for these irrationalities (because we use rational partial quo- 
tients, instead of integers), but between them, using the acceleration method of the 
previous section, we can find at the same time the approximations derived from the 
the secant method, the Newton method and the Halley method. 
We remember that a continued fraction is a representation of a real number a 
through a sequence of integers as follows: 

1 

a = a H , 

at H 

a-2 H : 

a 3 H 

where the integers ao,a\,... can be evaluated with the recurrence relations 
afc = [oik] 

1 . fc = 0,l,2,... , 

a/c+i = if ah is not an integer 

oik ~ a k 

for ao = a (see, e.g., [8J). A continued fraction can be expressed in a compact way 
using the notation [ao, ai, 0,2, a>3, ■ ■ ■]■ The finite continued fraction 

[ao, ...,a n ] = C n = — , n = 0, 1, 2, . . . , 

is a rational number and is called the n— th convergent of [<Xq, a\, 02, 03, • • .] and the 
aj's are called partial quotients. Furthermore, the sequences (p n )n=o an d (ln)n=o 
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are recusively defined by the equations 

)Pn 



O-nPn-l +Pn-2 
a n Qn-l + 



Vn > 2 , 



(4.1) 



and initial conditions po = ao, p\ = aoai + 1, qo = 1, q\ = a\ . In the following we 
will use rational numbers as partial quotients, instead of the usual integers. 

Remark 4.1. In general a continued fraction can have complex numbers as partial 
quotients. In this case given a real number there are no algorithms which provide 
the partial quotients. However, such a continued fraction can converge to a real 
number. In the following, we will study continued fractions of period 2 with rational 
partial quotients which are convergents to the root of larger modulus of a quadratic 
equation. For a deep study of the convergence of continued fractions, i.e., for an 
analytic theory of these objects, see [9] . 



Remark 4.2. A continued fraction with rational partial quotients 
an equivalent form as 



ao a± 
bo ' h ' 



has 



a Q 
to 



a i 



bib 2 



bob. 



a 2 



2 "3 



«3 



but the representation 



a 



ai 
h 



a 2 



1 



03 



is more suitable for the study of the convergents, as we will see soon. 



When we study the continued fraction 



ao ai a 2 
6 ' V b 2 ' "' 



, the equations (|4.ip become 



b„ 



-1 +Pn-2 
-1 + q n - 2 



Vn > 2 , (4.2) 



ao aoai 

with initial conditions po = 7—, p\ = — h 1, qo 

bo b h 



1, qx = — , providing two 



sequences of rationals . In the next Proposition we study how to determine the 
convergents through the ratio of two recurrent sequences of integers. 



Proposition 4.3. Given the continued fraction 



ao a\ 02 
bo' bi 6 2 ' "' 



, let (p n 



>n^0 and 



(q n )n=o be the sequences which provide the sequence of convergents (C n )^^ defined 
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recursively by equations If we consider the sequences (s n )^=0' {t n )n=Qj ( 



ii 



defined for all n > 2 

f s n = a n s n -\ + b n bn-is n ~2 
t n = a n t n -\ + b n b n -\t n -2 and initial conditions < 

— b n U n —\ 



nj n =0 



so = ao, si = aoai + bob\ 
to = 1, t\ = ai 
u = l 



then p r , 



b u r 



and q Tl 



t, 



for every n > 0, i.e., 



S„ 



b t r 



Vn > 0. 



Proof. We prove the thesis by induction. When n = 0, and n = 1 we have respec- 
tively 



Po 



ao so 



Pi 



bo bono 
aoa\ + bob\ s\ 



1 



to 



boh 



b U! 



qi = t- = — 



Now considering that 



Pn = J^Pn-l + Pn-2, n>2 
On 



by induction hypothesis we have 

a n S n -1 S n -2 _ a n S n -\U n -2 + b n U n -\S n -2 

b n bou n -i bou n -2 bob 

n'Un—l'an—2 
a n S n -lUn-2 + b n b n -iU n - 2 Sn—2 Sn 



b Q u n u n . 



boU r 



Similarly 



□ 



We finally focus our attention on a particular continued fraction of period 2. 
Precisely, we consider the continued fraction 



b b 

a c 



(4.3) 



for a,b,c G Z not zero. This continued fraction converges to the quadratic irra- 
tionality, of larger modulus, root of ax 2 — bx — c. Indeed, it is easy to check that 

b 1 b ca 2 2 2 

a = — | = — h 4=> aba + aca = b a + bc + aca 44> aa — ba — c = 0. 



a b 1 a ba + c 
-+ - 

c a 



From the previous Proposition we know that the convergents of fj4.3|) are determined 
by a linear recurrent sequence. Indeed, in this case the sequences s = (s„)^ and 
(t n )„^o of the Proposition 14.31 correspond to 

s = W{b, b 2 + ac, b, -ac), t = W(l, 6, 6, -oc) . 
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Thus, defining the sequence 

a = W(0,1, b, -ac) , 
we obtain that the convergents (C n ) of our continued fraction are 

C n = f^±L, Vn > . 
aa n+ i 

Clearly, they are not the best approximations of the quadratic irrationality, but since 
the convergents are determined by the ratio of the sequence a, we can apply our 
method of acceleration to such a sequence. In particular, we will see that properly 
accelerating the convergents sequence, we can find the approximations derived by 
Newton, Halley and secant method, which compare at the same time between these 
approximations. 



Theorem 4.4. Let us consider the continued fraction 



b b 

a c 



, whose sequence of 



convergents is (C n )^^ , and the sequence a = W(0, 1, b, — ac). If a is the root of 
larger modulus of at 2 — bt — c, then 

(1) (CF n+2 -i)n=o> are the approximations of a by the secant method; 

(2) (C2 n -i)n=Q> are the Newton approximations of a; 

(3) (CVi-i)^^, are the Halley approximations of a. 

Proof. We use the auxiliary sequence x = (x n )n5) defined by 

x n = , Vn > 2. 

C„-l 

We can directly apply Theorem 13. 1\ obtaining that (xF n+2 +i)n=o are ^ ne approxi- 
mations through the secant method of the root of larger modulus of t 2 — bt — ac. 
Considering the sequence x defined by 

x n = — , Vn > 2, 
a 



we have 



Cg n _ x Xg n _ 2 + ac aXg n _ 1 Xg n _ 2 + c 



« a x 9n _ 2 + xg n ^ - b ax 9n _ 2 + ax 9n ^ - b 
i.e. (x 9n )n=Q are the approximations through the secant method of a. Thus, re- 
membering that aC n = x n+ 2 for every n > 0, we have that {CF n+2 -i)n= are the 
approximations through the secant method of a. 
Applying Theorem 13.21 to the sequence x we obtain 



X2n-1 



x n + ac 



2x n - b 

i.e., X2, X3, X5, xg, . . ., X2 n +i, • • • are the Newton approximations for the root of 
larger modulus of t 2 — bt — ac. Now it is easy to see that x~2,x~3,X5, ...,X2«+i are 
the Newton approximations of a for every n > 2, i.e., (CV-i)^^ are the Newton 
approximations of a. 

Finally, using Theorem 13.21 again, with similar observations, it is possible to prove 
that (C3»i_i)^^ are the Halley approximations of a. □ 
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